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Abstract 

We consider different types of processes obtained by composing Brownian motion B(t), 
fractional Brownian motion Bn(t) and Cauchy processes C(t) in different manners. 

We study also multidimensional iterated processes in R d , like, for example, (Bi(|C(i)|), Bd(\C(t)\)) 
and (Ci(\C(t)\), ...,Cd(\C(t)\)) , deriving the corresponding partial differential equations satis- 
fied by their joint distribution. 

We show that many important partial differential equations, like wave equation, equation of 
vibration of rods, higher-order heat equation, are satisfied by the laws of the iterated processes 
considered in the work. 

Similarly we prove that some processes like C(|iJi(|.B2(— |-Bn+i(i)|—)l)l) are governed by 
fractional diffusion equations. 

Key words: Iterated Brownian motion; Fractional partial differential equations; Riesz 
fractional derivative; Fractional Brownian motion; Cauchy process; Wave equation; Vibration 
of rods. 

AMS classification: 60K99; 35Q99. 

1 Introduction 

The iterated Brownian motion is obtained by the composition of two independent Brownian motions 
Bi and £?2, as follows: 

J(t)=Bi(|Ba(*)|), *>0. (1.1) 

Recently this kind of processes has been studied by many authors: see, for example, Burdzy (1994), 
Khoshnevisan and Lewis (1996), Allouba and Zheng (2001), Allouba (2002), Beghin and Orsingher 
(2009). As far as the applications are concerned, it has been observed that the iterated Brownian 
motion I is suitable to describe diffusions in cracks (De Blassie (2004)) and many other physical 
phenomena. In particular, fractures on a rectangular slab can be viewed at as trajectories of a 
Brownian motion (see Chudnovsky and Kunin (1987)). The flow of a gas on a fracture can be 
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represented by a Brownian motion moving on a Brownian sample path and therefore this model 
produces an iterated Brownian motion (see, for example, Khoshnevisan and Lewis (1999)). 
In Orsingher and Zhao (1999) it is stated that the density of the process 



fY) X i' 

e 27 e 2r 



Pl (x,t)=2 ds (1.2) 

Jo \2ns \/2nt 

satisfies the fourth-order equation 

dp 1 d A p 1 gP 

^ (x ' i)= 2^^ + 2 72^^ <5(x) ' (L3) 

for x G R and i > 0. 

In De Blassie (2004) a similar result is proved under the initial condition p(x, 0) = f(x). 
Moreover it has been shown in Orsingher and Bcghin (2004) that the density (|1.2[) is a solution 
to the fractional equation 

(with initial condition p(x, 0) = S(x)), which is a limiting case of the fractional telegraph equation 

°P , 3/2 x d ^P 2 d2 P f- ttti , . . n 

"57 + 2 A TT-r = c -z-n, xeR,c,\,t>0, 
at Of? ox z 

for A, c — > oo, in such a way that c 2 /A — > 1. 

Generalizations of this result have been obtained in many directions. First of all by considering 
the n-times iterated Brownian motion 

I n (t) = B 1 (\B 2 (...\B n+1 (t)\...)\), (1.5) 

involving n + 1 independent one-dimensional Brownian motions. 

It has been shown in Orsingher and Bcghin (2009) that the law of (|1.5[) satisfies the following 
fractional equation 

««>'>°' <-» 

Bi(|B(t)l) 

, * >o, (1.7) 

where -B and £?i, Bd arc mutually independent Brownian motions, it is proved in Orsingher and 
Beghin (2009) that the joint law pj = pj(xi,X2, ...,Xd,t) is a solution to the fractional equation 

^f = ^7iE^I> i j eH,j = i 1 ...,d,i>0. (1.8) 



with p{x, 0) = 5(x). 

For the vector process 



2 



The constants appearing in equations (|1.4[) . (|1.6[) and (|1.80 have been chosen in such a way that the 
Brownian motions involved in the construction of the iterated processes possess standard distribu- 
tion. 

If we consider the slightly more general fractional equation 

ni d „ 2 

V^ = A2 E^' Xj €R, j = l,...,d,t>0, (1.9) 

3=1 J 

then the joint density of the vector process reads 

p d I (x 1 ,x 2 ,...,x d ,t) = 2 ===== H—=dw. (1.10) 

Result (jl.lOp shows that the components of the d-dimensional vector process (Bi(\B(i)\), Bd(\B(t)\)) 
are no longer independent and the parameter A enters into the variance of the "time process" 
B(t),t > 0. 

For the law of the vector process I d (t) we can also write the equation which is analogous to 

2 



j=i fc=i 

with initial condition p(cci, Xd 7 0) = IIfc=i H x k)- 

As far as the <i-dimensional case is concerned, along with the definition (| 1 . T[) of the process 
I d with components being independent one-dimensional iterated Brownian motions, another choice 
is to define the process as Bu)(\B(t)\),t > 0, with Br^ being a d-dimcnsional Brownian motion 
independent of the Brownian motion B{t). For the process I d (t) Allouba and Zheng (2001) and De 
Blassie (2004) show that the function u(t,x) = E x (f(I d {t))) = E{f(I d (t))\I$ = x) solves the initial 
value problem 

^-u(t, x) = & 2 u(t, x) + —j=Af(x), u(0, x) = f{x), (1.12) 
at V7rt 




for t > 0, x G R d , and 'sufficiently good' functions f(x), x € R d (which are supposed to be bounded 
with bounded Hlder continuous second derivatives). 

Moreover, when replacing the outer process Bu\ with a continuous Markov process X and 
considering Z(t) = X(\B(t)\), equation (|1 . 12[) remains valid for u(t,x) = E x (f(Z(t))) if we replace 
the Laplacian with the generator L x of the continuous semigroup associated with the Markov process 
X, that is: 

d 1 

— u(t, x) = L 2 x u{t, x) + —=L x f(x), u{0, x) = f(x), 
ot Wirt 



f G D(L X ), where D(L X ) is the domain of the operator L x . Moreover u(t,x) = E x (f(Z(t))) solves 
the fractional Cauchy problem 

■~-\-u(t,x) = L x u(t,x), u(0,x) = f(x), x G K rf , t > 
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(see, Baeumer, Meerschaert and Nane (2009)). 

A further extension is given in Baeumer, Meerschaert and Nane (2009) (see also Nane (2008)). 
Let X (t) — x + Xq (t) , where Xq (i) is a Lvy process in M. d starting at zero, and let L x be the 
generator of the semigroup p(x, t) = E x (f(X{t))). Then for any / G D{L X ) and any n = 2,3, 4, ... 
it was shown that both Cauchy problems for the higher-order partial differential equation 

a n— 1 ,1-1 

~i = J2w(I) L -f + L >' u{x,0) = f(x), xeR d ,t>0 (1.13) 

i—1 \n) 



and for the fractional equation 

d" u 

mi 



L x u, u(x, 0) = f{x), x e R d , t > 0, (1.14) 



have the same unique solution. This is given by u(x,t) = E x (f(X(EV n {t)))) , where £?V«(f) is the 
inverse Lvy subordinator of order —, that is, E 1 ' 71 ^) = inf{s : D(s) > t}, where D(s) is a stable 
subordinator with index —. 

For the Markov process with the inverse Lvy subordinator of order ^ one gets the same governing 
equations as in the case of a Markov process with a Brownian time, that is, these two processes have 
the same one-dimensional distributions. 

In the light of the above equivalence of (|1.13[) and 1)1.14)1 (and inspired by the result of Orsingher 
and Beghin (2009), equation (|1.6|l ). the following PDE connection of n-iterated Brownian motion 
appears in Nane (2008): the unique solution of the Cauchy problems (|1.13|) and (|1.14|) with n = 2 k 
is given by u(x,t) = E x (f(X(Ik-i{t)))), where X is the Lvy process described above and I kit) is 
fc-times iterated Brownian motion (|1.5|) . 

For the law of the n-times iterated Brownian motion, we can write an analog of equations (|1.3|) 
and (TETT)) : 

i=E ^7T){o^) S ^ + {o^) p ' PM) = *(*), ^R- *>0, (1.15) 

where m = 2". 

For the fractional diffusion equation of order 1/n, iieN, 

dip 1 d 2 p 

^| = -^|, , e R, t >0, (1.16) 

with initial condition pi/ n (x, t) = 5(x), it has been shown in Beghin and Orsingher (2003) that the 
solution coincides with the density of the process 

h/n(t) = B m Gj(t) j , n > 1, (1.17) 

where (Gi(i), G n (t)) is a random variable with joint distribution 



p(wi,...,w n -i) = i r e " w 2 -w™_l, Wj > 0, 1 < j < n- 1. (1.18) 

(2tT) 2 y/t 
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Taking into account the result of Baeumer, Meerschaert and Nane (2009) on the solution to 
the equation (|1.14p . we can conclude that the process (|1.17[) and the process B (£ 1/n (i)) , that is a 
Brownian motion subordinated to the inverse Lvy process of order i, have the same one-dimensional 
distributions, which can be written down in a closed form. 

In D'Ovidio and Orsingher (2009) various types of processes obtained by composing independent 
fractional Brownian motions Bn 1 (t) and Bn 2 (t),t > (with Hurst parameters respectively equal to 
Hi,H 2 £ (0, 1)) have been studied. In particular, for the density phxH-2 (#j *) of Bh x (\Bh 2 (t)\) it has 
been shown that the governing equation has the following structure: 

(1 + + t 2 gf = H\Hl fjtef? + x 2 Sl , X e R, t > 0. 



dt dt 2 1 1 \ dx dx , 

In this paper we extend some of these results in several directions, by considering different compo- 
sitions of Brownian motions, fractional Brownian motions and Cauchy processes. 

We start by studying the iterated Brownian motion defined in in the case where B\ is 

endowed by drift, giving the two governing equations, the fractional and fourth-order one (sec (|2.4[) 
below). 

Another direction of our research concerns processes obtained by combining Brownian motions 
(possibly fracional Brownian motions) and Cauchy processes. In particular wc study the multidi- 
mensional vector processes [B x (| C(i) | ),..., B d (| C{t) \ )) and {C x (| B(t) | ),..., C d (| B(t) | )) , which involve 
independent Brownian motions Bj(t),t > and Cauchy processes Cj(t) 7 t > 0. 

For d = 1 it is has been proved in D'Ovidio and Orsingher (2009) that the law of B(\C(t)\), t > 0, 
given by 

2 f°° e~i t 

PBc(x,t) = - -== 7 ■ 7 ds, x£R,t>0 (1.19) 



7T 

satisfies the following fourth-order equation 

d 2 v 1 d A v 1 d 2 

with initial condition p(x,0) = S(x). We extend this result to the case where the Brownian motion 
is substituted by the fractional Brownian motion BH(t),t > 0, with Hurst parameter H £ (0,1). 
Wc show that the density function of Bjj(\C(t)\),t > resolves the following equation 



dt 2 



2Ht * SI M 



with p(x, 0) = 6(x). 

Also the case of a Cauchy process with a random time represented by an iterated Brownian 
motion is taken into account and the corresponding fractional governing equation is derived. 

The iterated Cauchy process, defined as Jcc{t) = Ci(|C2(i)|), t > has been considered in 
D'Ovidio and Orsingher (2009) and its connection with the wave equation 

has been established. 
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We will show here that when the process C± is endowed with a position parameter a / 0, the 
distribution of the iterated Cauchy process satisfies the following wave equation 

d 2 v d 2 v 1 

W = ld-^-ar> *€R,t>0. (1.22) 

Finally we consider the <i-dimensional vector (Ci(\C(t)\), Cd(\C(t)\)) and obtain the equation 
governed by its density, i.e. 

2 r+oo d ^ ^ 
P d cc (x 1 ,...,X d ,t) = ^1 J o II s 2^ T 2 t 2 + s 2 dS ' ( L23 ) 

where the use of partial Riesz fractional derivatives is required. 

2 Iterated Brownian motion with drift 

If the processes composing the iterated Brownian motions possess drift, the fractional equations 
and the higher-order equations governing the distributions are somewhat different. We start by 
considering the process I^(t) = B^(\B 2 (t)\),t > 0, with law 

e 2* 



p1(xA) = 2 / == =ds. (2.1) 

7 Jo v^s V^rf 

It has been shown in Bcghin and Orsinghcr (2009) that (|2.1[) solves the following fractional 
equation of order v = 1/2 : 

dip 1 d 2 p u dp _ .„ „. 

e$=w*i£--%£ ieE - t>0 - (2 - 2) 

As a check we evaluate the Laplace- Fourier transform of (|2.1I) : 



J-00 \ JO 



e 2t 



+00 / 

e W* 2 / c " ^-==ds ) dxdt (2.3) 



27TS V27rf 



00 . , p-^v^ 



= 2 e~* p S+W1S — ^ds 







The previous result coincides with the Fourier-Laplace transform of 
We show in the following theorem that, as in the case where fi = 0, the density of I^(t) satisfies 
also a fourth-order equation. 

Theorem 2.1 The density of the process I^(t) = B^(\B2(t)\),t > given in \2. 1}) is a solution to 
the following equation 

B , , 1 / 1 d 2 d \ 2 , , 1 fid 2 d 



dt PM = 2[2dx- 2 ^dx PM + 7m\2d^ 2 ~^ 5 ^ M 
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with initial condition p(x,0) = 5(x). 

Proof By taking the first-order derivative of (|2.ip . we have that 



e 2 S (j e 
\/2ns dt J2ni 



ds 



d 2 e" 



2tt.s 9s 2 ^/27^^ 
[integrating by parts] 



9 e 



d e 



o 3s v / 27rs ds y/2ni 



_s£ 

o e 2s e 2t 



9 s 
1 



27TS V27rf 



f/.S 



oo a2 



9 2 e" 



<9s 2 



^TTS V27Tt 



ids 



1 d 2 



27rt V 2 dx 2 ^'dx 



5(x) 



d fl d 2 



1 /l 9 2 







dx 



1 



_9_ 

ds V 2 <9x 2 ^ dx 



1 d 2 



2 Ua?"^ I W' t)+ ^loX2-^z)*( 



27rf V 2 dx 2 



dx 



e 2* 



^tts V2irt 



ds 



(2.5) 



□ 



Remark 2.1 The differential operator appearing in the fourth-order equation (|2.4p is the formal 
square of the operator appearing in the fractional equation (|2.2p . In the special case where [i = 0, 
we obtain again equation (jl.3p . 



We consider now the case where the process representing the "time" possesses drift, then the law 
of the iterated process = > 0, reads 



qj(x,t) 



C{t) Jo 



e 2 S e 



<2lZS V27Tt 



-ds, 



where 



C(t)=Pr{S^(t)>0} 



e 2 



dy 



is the normalizing factor for the density of \B% (t)|. 
The Fourier transform of qj(x,t) becomes 



(2.6) 



qj(x, t)dx 



-il 
e 2 



(y/tw + flt)- 



C{t) J_^ t 



2tt 



-dw. 



(2.7) 



The evaluation of the Laplace transform of (|2.7p poses serious problems. For this reason the case 
where the process representing the "time" possesses drift is not further developed here. 
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3 Iterated processes involving the Cauchy process 

The transition function of a centered Cauchy process C(t),t > is given by 

Pc{x, t) = ^2^2) » t > (3.1) 

and it is well-known that pc is a solution to the Laplace equation 

[d^ + ^) PcM =°- (3 ' 2) 
Furthermore pc is a solution to the following space-fractional equation 

d d 

—p c (x,t) = -——p c (x,t), (3.3) 
ot o\x\ 



where 



d\x\ J[X) ndx 



- m dy _rm dy 



(3.4) 



'-co x-y j x y-x 

is a special case (for v = 1) of the Riesz fractional space-derivative and possesses Fourier transform 
equal to 

f + OO O /-+0O 

e i/5a -5i-r/(*)<fc = 1^1 / ^ X f(*)dx = (3.5) 

By composing the standard Brownian motion and the Cauchy process, we obtain the following 
new processes: 

J BC (t) = B(\C(t)\), t>0 (3.6) 

and 

J C B(t)=C(\B(t)\), t>0, (3.7) 

where C and i? are mutually independent. 

An alternative definition of the first-order fractional derivative is given in Saichev and Zaslavsky 
(1997) (see formula (A. 39)) and reads 

d , 1 f + °° f(x-y)-2f(x) + f(x + y) j 

M f{x)= -,L ? dy - m 

This is a special case of 

d v , _ 1 /• + °° /(a; - y) - 2/(x) + /(i + y) 



i\x\^ {x) 2r(-^)cos^i 



2 



dy, (3.9) 

for < v < 1. We note that (|3.9[) reduces to (|3.8[) for u = 1, because, by the reflection formula of 
the Gamma function, the constant can be rewritten as 

i i r(i + z/) 



2r(-zy)cos^f 2 cos 



7Ti/ 7T 



2 sin-TTi^ 

r(l + i/) sin^ 
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which, for v = 1, yields —1/tt. 

We can convert definition Q3.4p into (|3 .8[) by means of integrations by parts, after suitable changes 
of variables. 

As a further check we can show that the Fourier transform of p. 81) coincides with (|3.5[) : 



+00 



+00 



f(x - y) - 2f(x) + f(x + y) 



dy dx 



f{w)e lPw dw 



— OO 



y 



+00 



+00 



f(w)e l/3w dw 

'0 y 



dy+- / f(x)e ipx dx 



-dy 



V 



;dy 



dy. 



(3.10) 



The last integral can be evaluated as follows: 



+00 e i/3y _ j /.+00 g-i^a _ ^ 



V 

+00 



-dy 



V Jo 

" l/3w dw 



+00 



-dy 



-^dy / i(3e if3w dw- / -jdy / ipe'^dw 



r Jo 



i/3 / e^dtu / —dy-i/3 / e"^ 10 ^ / — dy 

Jw y Jo Jw y 



2d 



+00 ij3w 



-2/3 



+00 



2wi 
sin /3ui 



dw-2/3 



+00 e -ip w 



2wi 



-dm 



dw = — 7T I 



which, inserted into (|3.10[) . gives p.5j) . 

It has been shown in D'Ovidio and Orsinghcr (2009) that the transition density of the process 
JBc{t) = B(\C(t)\), t > 0, which is given by 



, , 2 f 00 e"^ t 

PBC{X,t) = - — == 

7T Jo V 27TS ^ + 



ds, 2; G R, i > 



(3.11) 



satisfies the following equation 



<9 2 p 1 <9 4 p 1 d 2 r , , 



with initial condition p(x,0) = 5(x). 

Remark 3.1 In mathematical physics the equation 



(3.12) 



= —K- 
dt 2 dx* 



(3.13) 
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represents the vibration of rods (see Elmore and Heald (1969), p. 116). Equation (|3 . 1 3[) coincides 
with ff3TT2|) for x ^ 0. 



Remark 3.2 The density (|3.11[) can be worked out in an equivalent form, by applying the subor- 
dinating relationship 



te 2» 



-.dw. 



ir(t 2 + s 2 ) J V2ttw 3 
Thus the transition density of the process lBc{t) = B(\C(t)\), t > can be rewritten as 



PBc{x,t) = 2 



e 2s 



e a™ te 2 « 



27rs Wo \J2ttw a/27TW 3 



zdw ds. 



(3.14) 



Formula (|3.14j) corresponds to the law of an iterated Brownian motion I(t) = Bi(\B 2 {t)\), taken 
at a random time T(t), which coincides in distribution with the first-passage time of a standard 
Brownian motion. Thus (|3.14[) coincides with the density of the process I(T(t)),t > 0, where 
T(t) = inf(s : B(s) = t). 

We consider now the composition of a fractional Brownian motion with a Cauchy process, 
J BH c(t) = B H (\C(t)\),t>0. 

Theorem 3.1 The density of JB H c{t)it > 0? which is given by 



PB H c(x,t) = - 

TT 



ds, 



satisfies the following equation 
t 



dt 2 



d_ 

dx" 



8 2 

2 U 2 



H(H-l)— x-H-— x 

ox 2 



2Ht d 2 

p -—dx^ S{x)1 {°<»<iY 



with initial condition p(x, 0) = S(x). 

Proof We take the second order time-derivative of (|3.15[) . that is 



(3.15) 



(3.16) 



d 2 

gpPB H c(x,t) 



2 f +QO e"^ d 2 

TT 



II 



t 



V r 2r7s m dt 2 \t 2 + s 2 
t 



ds 



2 

TT 

2d_ 

K~ds~ \ V / 2^2H I t 2 + s 2 



^rTs^ds 2 \t 2 + s 2 

s=-\-oo 

t 



ds 



s=0 



ds 2 \ V2^2H I t 2 + s 2 



ds. 



(3.17) 
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The first term in (|3.17[) can be worked out by considering that the transition density of the fractional 
Brownian motion pn = pn{x,t) satisfies the following equation 

#P = Ht 2H-i d 2 P 
dt dx 2 



and it becomes 



2H ts 



2H-1 Q2 



e 2» 



IT t 2 + S 2 dx 2 \ ^/2^H 



s=0 



\<H<1 



The second term in (|3.17p can be evaluated as follows 



+°° d_ 
~d~s 



Hs 2 "- 1 & 



e 2 S 



3a; 2 \ y / 2irs 2H t 2 + s 



Q<H<\ 



■ds 



(3.18) 



(3.19) 



2H(2H - 1) 



n 2H-2 u 



t 



2H 2 



dx 2 \ yph^x J t 2 + S 2 
t 



ds + 



e 



dx 4 \ V^H" / t 2 + S 



-ds 



2H(2H - 1) d 
tt dx 



+oo 



e 2s 2 -£ 



t 



2H 2 



2H-2 u 



e 2 



V2^ 2W t 2 + s 2 
t 



ds 



dx 3 \ V2^s 2W / t 2 + s : 



-ds. 



We rewrite the last line in f|3 . 19f) as 

2H 2 



+oo o2 
g2H-2 ° 



dx 2 \ V2^s 2W / t 2 + S 



;ds + 



(3.20) 



2H 2 



2 /•+oo 



t 



2H 2 



2 r+°° 



dx 2 
,d_ 

dx \ " y/2ns 2H J t 2 + s : 



V2^ 2W ) t 2 + s 2 
f 



ds 



■ds 



2H 2 



2 r+°° 



dx 2 



V2^ 2ir / t 2 + s 



:ds. 



By inserting (|3.20[) into f|3 . 19[) the second term of (|3.17|) takes the form 



2H(H-l) d 2H 2 d 2 2 

X — — ttX 



dx 



7r dx 2 



V2^s 2lT t 2 + s' 



■ds. 



(3.21) 
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The integral appearing in (|3.2ip can be calculated by observing that 

2 

r+oo — - 



t 2 + s 2 - s 2 



In , , 1 
= --p2VB H c(x,t) + ~ 

The integral in p.22[) can be explicitly calculated as 



72^1/ t 2 + S 2 



ds 



V^S 2 " 



ds. 



+00 



V^s 2 " 



ds 



.r 



l/H 



2 1/2H 

e 



-1/2H 



-y x -2/H y l/H X l/H 



1 

2H 



22«^ 



1 



1 



xtt +1 H^ \2H 2 / 
which, inserted in p.22[) , permits us to rewrite Q3.2ip as follows: 



2H(H - 1) d 2H 2 d 2 



7T dx 

H(H-l) 

r 



dx 2 



1 7T 



2 27T- 



1 

22TT(H - 1) 



H(H-l) 
2w(g-l) r 



1 

72 



ff(fT-l) 



dx 


H 2 


—x 2 

dx 2 


(±. 






\2H 


*\ 




d 
dx 


H 2 


dx 2 


f_L. 




j X~TT~ 


\2i? 


f * 


a 

dx 


H 2 


dx 2 



PB H C + 

1 \ 2m if 

« — r 

/ ir^/nt 

PB H C + 

1 2th(H-1) 



1 
2i7 



(P_ 

dx 2 



-+1 



1 
277 



a; ff 



By putting together (|3.23l) and (|3.18[) we obtain equation (13.161) . 



(3.22) 



(3.23) 



□ 



Remark 3.3 In order to check that for H = i equation (|3.16[) can be converted into (|3.12l) . it is 
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necessary to perform the following calculations: 

t 



J*. 

dx 4 



+°° e - ^ 



o 



^t 2 



ds 



dx 3 
2 d 2 

7T dx 2 



2 

— X 

7T 



e 2 3 



/2^< 2 



1 e 2 S 



^i 2 



(3.24) 



In analogy with (|3.22[) , we can write 



Furthermore 



1 e 2 S 



1 K / s 1 

= —p-jPBCKXyt) + - 



+oo 



1 e 2s 



'27TS 



(3.25) 



9a; 2 \ 

dx 

_d_ 

dx 



e 2s 



sVZns t 2 + s- 



-</.s 



/2^t 2 



-ds 



1 7T . . 1 

~ — -PBC{X,t) + - 



t 2 2 J 



+ 00 



1 e" 



/27TS 



By inserting ([335) and (|3T26|) into (pT24|) we have that 
S 4 



dx 4 



PBc{x,t) 



Id 2 Id 

J2Q^2 (xpacfat)) - (x PBC (x,t)) + 



2_d*_ 

irt dx 2 



'2ns 



2 d / 
=ds H ~ — x 



7rf eta 



27TS 



zds 



1 d 2 . Id 

-pQ^ (x 2 PB C {x,t)) - -— (xp BC (x,t)) 



2 



( 9 2 



<9a; 2 dx 



e 2s 



/27TS 



zds. 



(3.26) 



The calculations carried out in Theorem 3.1 permit us to conclude that, for x =/= 0, 
d 4 i g 2 19 d 2 

g^PBc(x,t) = -p-g^ {x 2 p BC (x,t)) - - — {xp BC {x,t)) = ~—p BC (x,t). 

Therefore, despite the formal diversity of (|3.12[) and (|3.16[) , the law of B(\C(t)\) satisfies also equation 
©T6]). 



13 



We consider now the <i-dimensional generalization of the results given in p. lip and (|3.12p : let 
us define the vector process 



J d BC (t) 



Bx(\C(t)\) 



B d (\C(t)\) 



t > 0, 



where B±, ...,Bd are Brownian motions independent from each other and from the Cauchy process 
C(t),t > 0. Then the following result holds. 

Theorem 3.2 The joint density of Jg C {t),t > 0, which is given by 



+ 00 



n 



e 2s t 
/2ns t 2 + s 2 



ds, 



(3.27) 



is a solution to the following equation 



\J=1 3 ) 3=1 3 



(3.28) 



with initial condition p(xi, x<j, 0) = Ylj=i^( x j)- 

Proof By taking the second order time-derivative of (|3.27[> we get, by means of two integrations by 
parts, that 



92 d 

~q£2PbC\ x 1> ■■■i x d,t) 



ds 2 



t 



/2ns j t 2 + s 2 



ds + 



--in 

nds I 11 



t 



J =1 



/2tts I t 2 + s 2 



+ 00 



(3.29) 



2 [ 1 y>_^_ 



+00 



n 

3=1 



e 2s t 
/2^t 2 + s 2 



ds + 



y— in 

nt 2 + s 2 Z-< dx 2 11 

3=1 3 Vj=i 



which coincides with ([3.280 . 



□ 
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Remark 3.4 In view of Remark 3.2, we note that the vector process Jgc(t) is equivalent in distri- 
bution to the following d-dimensional process: 



h{T{t)) 



h(T(t)) 



t > 0, 



where Ij, j = 1, ...7i are independent iterated Brownian motions and T(t) coincides with the first- 
passage time of a standard Brownian motion. 



We focus now our attention on the process 

J CB {t) = C{\B{t)\), 
and prove that its density 

f s 2 f +oc s e~£ , 

satisfies a non-homogeneous backward heat-equation. 
Theorem 3.3 The density \3.30\) is a solution to the following equation 

1 



t > 0, 



dp 
di 



1 d 2 P 

2dx* + ^V^Ft 1 



x € R, t > 0, 



with initial condition p(x,0) = S(x). 

Proof By taking the time-derivative of pcb = PcB(x,t), we have that 



dp 

~dt 



8 



x 2 dt ' 

d 2 



/2TTt J 

2 

e~ 2* 



ds 



s 2 + x 2 ds 2 \ 

+°° d_ f s_ 

ds \ s 2 + x 
ds 2 



ds 



I^TTt 



-ds. 



which coincides with (13.311) . 



(3.30) 



(3.31) 



□ 



The previous results can be further extended to the case where the "time process" is represented 
by a n-times iterated, instead of standard, Brownian motion. We prove that the density of the 
process 

C(|I n (t)|) = C(|B 1 (|B 2 (...|B fH . 1 (f)|...)|)|), * > 
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is a solution to a non-homogeneous fractional equation. 

Theorem 3.4 The density of the process C(\I n (t)\),t > 0, which is given by 



2 r+oo 

Pci(x,t) = - 

n Jo 



;P„(s,t)ds, 



where 



p n (x,t) = 2 r 



-oo r-\-oo 



e 2 ™i e 2 ™2 e 2* 

dwi y= — =dw2--- ^ — - dw ni 



satisfies the following equation 

d^p „ 1 9 d 2 p 2 



r(-|) 



9x 2 



- 3 / 2 r(-^ 



.r- 



'27rt 



iel, i > 0, 



wzi/i initial condition p(x,0) = S(x). 

Proof We start by taking the fractional time-derivative of order 1/2™ of the density (|3.32[) 



d 2" 



-p C i(x,t) 



8t& 

2 r + °° s d 

7T Jo S 2 + X 2 dt^ 
2iW-l /.+00 s Q2 



7T JO 

2^-1 



—p n (s,t)ds 



Pn{s,t)ds 



,■2 



ds 2 



8 



7T S 2 + X 2 8s 

25W-1 ^+00 5 



Pn(s,t) 



jo ds \s 2 + x 2 J ds 



8 



p n (s,t)ds 



22-" 1 5 

7r 9s \s 2 + a; 2 

2^-1 ,+00 g2 



Pn(s,t) 



fx 



+ 



7T Jg 8s 2 \S 2 + X 



Pn(s,t)ds 



22W- 1 1 



1 2 a 2 



7T X 



r p„(0,t)-2^- 2 ^-2 Pc/ (x,t) 
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We concentrate now on the first term and evaluate 



Pn(0,t) = 2" / _ dWi / _ dw 2 - / nr= dw n 

J ^2nwi J V 2 ™2 Jo V27rw„v27rt 



2 n -ir(^) f + °° _i »| / ,+0 ° 



,4; 



w 2 4 e 2 "s di«2... / 



e 2t 



1 1 . 



2tt 



1\ i i /3\ i i fl 1 



2-3-ir - 2-3-sr - ...2-2-2-r / — = dw 



n+1 W V 8 /'" V 2 2n J Jo \ft 



2™ ii /l\ i i /3\ i i (I 1 \ _i _j / 1 1 

_2^2-4r - 2^5-sr - ...2"'-^r 2 



n+1 \4/ \8/"' V 2 2" 7 V2 2" +1 



2^)' 



2 -l+iSTT /J 1\ /l 1 \ /l 1 \ / 1 1 



tt^t 1 \2 2 2 / V 2 2 3 / "' V 2 2 " / V2 2 ri+1 

[by the duplication property of the Gamma function] 



2 ^SnTT^^^^^^ r( 



■1+1 

7T 2 



r(-^r) 



1 V 2, 



r(-^ 



which, multiplied by the constant appearing in (|3.34[) . gives the final form of equation p.33[) . □ 
In the d-dimensional case we consider the vector process 

" Ci(\B(t)\) 

t>0 (3.35) 

C d (\B(t)\) 

and obtain the governing fractional equation in the following theorem. 
Theorem 3.5 The joint probability law of the process defined in US. 35\) reads 

2 r+°° JL s e — §| 

PcB ( Xl ,..., Xd ,t) = - d j o n^—^ * ( 3 - 36 ) 

j — 1 j 

and satisfies, for al > 1 7 the following fractional equation 

d 

with initial condition p(x\, Xd, 0) = Yi ^( x j)- 
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Proof The time-derivative of (|3.36[) can be evaluated as follows: 
9 



rj (>-\-00 « 

7^ l n 



die at \ 

as 



L s 2 + a; 2 9* I V2?rf y 



i/ n 



s d 2 I e at 



- 2 + x 2 



" ./o s ^ + ^ s2 V V27rt 



ds 



--in 



e 2* 



7T d 9s \ 1J - S 2 + X 2 J V27rf 



*d /. 9s 2 ( jl s 2 + x 2 | VSrf 



f/.s 



e 2* 



1 d d s d ( s 

TT d ^ n S 2 + X 2 8.8 \S 2 +X 2 J y/2^l 
k=l j = l J \ J / 



+oo q2 



(3.38) 



n 



n d J ds 2 \ 11 s 2 + x 2 I ^J 7 2^i 

\3~ / 



ds. 



The first term in the last member is equal to zero (unlike the one-dimensional case) and this makes 
equation (|3.37p homogeneous. 
Since 



d d 

En 







ds \ 1 1 s 2 + x 2 I 1 1 s 2 + x 2 ds \s 2 + xlj 



the second-order derivative becomes 



Q2 I d 



ds 2 



n 



s 2 + x 2 



E 

fc=i 



E 



d 



n 



9s \ s 2 + 7 9s \ s 2 + a; 



,2 i -2 



■? / 1=1 S +X l 



+ 



n 



ds 2 \s 2 + x 2 J j^s 2 + 



(3.39) 
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In view of p.2[) and p.3[) , we can rewrite (|3.39p as 



ds 2 1 n 



s 2 + x 2 



= E 



k=l 



5 



n 



J= t ^l^fcl V s2 + a; fe/ ^kjl \ s2 + a; i / 



l^k 

d d „ „ d 

O O -n- S 



(3.40) 



EE aK | 5 |^|ll s 2 +a; 2 



n 



e dx 2 ( n 



fe=i j 



k=l 



k \l- 



1=1 



1 i 2 
S z + xf 



By inserting (|3.40[) in (|3.38p . we arrive at equation (|3.37[) . by applying the commutative property 
of the fractional derivative (13.81). Therefore we show that 



d d 



d d 



d 2 



d\xj\d\x k \ d\x k \d\xj\ d\x k \d\xj[ 
The second-order fractional derivative reads 
d 2 

d\yW\ f{x,y) 

_d_(l f +Q ° f(x-t,y)-2f(x,y) + f(x + t,y) 

d\y\UJ t 2 



(3.41) 



dt 



+ °° dz f + °° f(x -t,y-z)- 2f(x -t,y) + f(x -t,y + z) ^ 
" Jo z2 t 2 

-— r°°dz r°° f{x > y ~ z) ~ 2f{x > y) + f{x > y+z) dt 



z 2 t 2 

-— [ +oc 'dz [ + °° f{x tlii^itz 2 ^ x +M+ ^ x + t > y + z \ t 



z 2 t 2 

It is easy to check that giJLui f( x ; y) produces the same result and thus the commutativity of the 
second-order fractional derivative holds. □ 



We consider now the last case of composition of Cauchy processes: it can be seen that the density 
of the process 

C?(\C 2 (t)\) t > 0, 

where the external process is endowed with a position parameter a 6 1, is a solution to a non- 
homogeneous wave equation. Indeed, by suitably adapting the proof of theorem 4.1 in D'Ovidio and 
Orsingher (2009), it is easy to check that 



Pcc(x,t) 



t 



s 2 + (x- a) 2 t 2 + s 2 



ds 
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is a solution to 

In the d-dimensional case the iterated Cauchy process can be defined as follows 

Ci(\C(t)\) 

4 c {t) = { ... , t>0 (3.43) 
C d (\C(t)\) 

where C\...Cd and C are mutually independent, standard Cauchy processes. 
Theorem 3.6 The density of J^ c (t),t > 0, which can be expressed as 

2 r+oo d s 1 

p d cc ( Xl ,...,X dl t) = - m II s 2 +x 2 t 2 +s 2 ds ( 3 - 44 ) 

j — l 3 

and satisfies, for al > 1 7 the following fractional equation 

w = 5 - 2 1 ? otSky Xj £R > j = 11 d ' ' > °' (3 - 45) 

A, — X A. — 2 j — 1 

d 

with initial condition p{x\, Xd, 0) = J} <5(xj). 

j=i 

Proof The second-time derivative of (|3.44j) can be evaluated by adapting the proof of the previous 
theorem, as follows 

92 d , \ 

2 /- +00 A g d 2 / t \ 

7T d+1 Jo S 2 + X 2 <9i 2 V 2 + S 2 / S 



n s 2 +x 2 9s 2 U 2 +s 2 J ds 



s 



n 



7r d+1 t 2 + s 2 ds \ i - i -s 2 + x 2 



2 /"+ 00 t d 2 



d 

S 



(I 



= [byd339l)] 



E fef , .... x* t) - 2 E E fli^ia, | P&o(»i. -.***)• 

A* — 1 At — 2 j — 1 
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□ 



Remark 3.5 The density of Jcc(*)j t > 

can be expressed in the following alternative form 
d 2t t 

dote, .... x d , t) = n -r^^) ln N ' (3 ' 46) 

as can be inferred from the calculations leading to theorem 4.1 of D'Ovidio and Orsingher (2009). 
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